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Abstract 

The uncertainties in estimating the hadronic production of the B c meson are studied under 
the framework of the complete approach of the perturbative QCD and the gluon-gluon fusion 
mechanism. Quantitative comparisons on the production at TEVATRON and LHC are made. 
Considering the detectors at TEVATRON and LHC, we have also estimated the production with 
proper kinematic cuts. Based on the results, we conclude that the experimental studies of the B c 
meson at the two colliders will be complimentary and stimulative. We find that as cm. energy is in- 
creasing from RUN-I to RUN-II at TEVATRON, the production cross-section increases about 20%. 
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I. INTRODUCTION 



-B c -physics is attracting more and wide interests due to experimental 



and theoretical ones 
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(especially the experimental discovery |3|]). Since the production cross section is relatively 
small in comparison with the production of the other heavy mesons (B, D etc) and heavy 



quarkonia (J/ip, T etc), only at high energy hadronic colliders with 



one collect enough B c events for experimental studies 
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111 ]. Considering 



the potential applications to the experimental feasibility studies on the topic, a computer 



program for hadronic generating the B c meson, BCVEGPY, has been completed in Ref . 23] , 
which has been written in Fortran and can be conveniently implemented into PYTHIA 24] . 

The hadronic production of the _B c -states with different spin has been estimated by a 
lot of authors, and all are mainly based on the gluon-gluon fusion mechanism, i.e. taking 
the hard subprocess gg^B c + b + cas dominant one. The authors _of 
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111 did it in the 



estimate in terms of the 'fragmentation approach', and the authors of 

so-called 'complete approach'. The results of the two approaches agree essentially within 

I I 

theoretical uncertainties, especially, at high B c transverse momentum P? jlOj and when 
the fragmentation approach additionally takes into account the contributions from gluon 
fragmentation Q]. The fragmentation approach is appropriate if one is only interested in 
the information of the produced B c itself, and its accuracy is improving with the increasing 
of Pt [ill- Comparatively, the advantage of the complete approach is that it can also retain 
the information of the accompanied c and b quarks (jets). For experiments, to retain more 
information of the events is more relevant. Therefore, in BCVEGPY |23j], the complete 
calculation approach has been adopted. 

Moreover according to perturbative QCD (pQCD), in addition to the gluon-gluon fu- 
sion mechanism, there are several different mechanisms for the production, such as that via 
the quark-antiquark annihilation subprocess qq —>■ B c + b + c, and that via the color octet 
mechanism and etc.. In Ref. sj], it is shown that the contributions to the production from 
quark-antiquark annihilation are much smaller than those from gluon-gluon fusion. It is 
mainly because of the fact that the relevant luminosity for gluon-gluon fusion is a product 
of the parton distribution functions (PDFs) for gluon components in the colliding hadrons, 
whereas 6that for quark-antiquark annihilation is a product of the quark and anti-quark 
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distribution functions respectively while in a high energy collider, such as TEVATRON (pp 
collision VS = 2 TeV) at Fermilab and LHC (pp collision = 14 TeV) at CERN, the 
luminosity for gluon-gluon fusion due to parton distribution functions (PDFsJ at the most 
kinematic region is much higher than that for quark- ant iquark annihilation [23] - Further- 
more, in the annihilation subprocess, there is an additional S-channel suppression from the 
virtual gluon propagator, which will also make contributions to the production cross section 
small. In the subsection II-C we will make a brief comparison between those two mecha- 
nisms. According to non- relativistic QCD (NRQCD) 25], for a double heavy meson, such as 
J/i/j,r] c , B c , ■ ■ •, there may be additional color-octet components which may contribute the 
production substantially. The magnitude of the color-octet components may be estimated 
with the NRQCD velocity scaling rule (v ~ a s (mB c ) ~ 0.1). One may find that there is 
no enhancement in the hadronic production (in contract to of the hidden flavored double 
heavy objects such as J/ip, ip'the production which can gain one order or more in a s via 
the color octet components) and the color octet components in B c are definitely smaller 
than the color singlet ones (with the same input parameters, the total contributions to the 
hadronic production is ^ v ~ v 2 of the color-singlet one), therefore, the contributions from 
color octet components of B c are ignorable at leading order. 

The contributions to the production from the color octet and those from the quark- 
antiquark annihilation mechanism were neglected in Refs-P. [loL 111 . 12, 13]. In BCVEGPY 
[23d and thus in the present paper, the mechanisms via color-octet for the B c production are 
ignored. 



Although the estimated values for the hadronic production of B c in literature 
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121 ] are consist within the theoretical uncertainties, in a few cases the differences on the 
production can reach up to almost one magnitude order. Moreover, RUN-II of TEVATRON 
is taking data, LHC is under constructing, and various experimental feasibility studies of 
B c are in progress, thus we think to know the theoretical uncertainties quantitatively in 
estimates of the B c production, in addition, the precise comparisons of the production at 
LHC and TEVATRON are interesting. 

The newly developed generator BCVEGPY has been well tested by carefully comparing 
the results with those in earlier references |8(. It, having been implemented into PYTHIA, 
can be conveniently applied to simulate B c events (Monte Carlo simulation). In the pa- 
per, we highlight the uncertainties of the dominant gluon-gluon fusion mechanism, make a 



3 



brief comparison on the gluon-gluon fusion mechanism and the quark- ant iquark annihilation 
mechanism but precise comparisons on the production at TEVATRON and LHC. Although 
some certain uncertainties will be able to understand (control) better when the next to 
leading order (NLO) pQCD calculation is achieved, considering that the NLO calculation 
cannot be available soon due to its complicatedness, we restrict ourselves here to examine 
the uncertainties only 'up to' the lowest order. The uncertainties that we will examine in 
the paper, include the variations about a s -running, the choices of the factorization energy 
scale, the various versions of parton distribution functions (PDFs), the values of the bound 
state parameters and etc.. To be useful references, the production with possible kinematic 
cuts, which roughly 'match' the detectors at TEVATRON and at LHC, is also investigated. 

It is known that at high energy hadronic colliders, numerous B c events may be produced, 
moreover due to the fact that the B r m eson has a quite large branching ratio to decay into a 
B s meson (several tens percent) 16. \l\ \\\. thus, in principle, copious and precisely tagged 
B s mesons at production position may be collected through the inclusive decays of B c at the 
hadronic colliders (B c (B c ) is the charged meson, as long as the charge of B c (B c ) and/or 
its decay products is determined, the produced B s (B s ) mesons through the decays of B c 
(B c ) are tagged precisely at the decay vertex of B c (B c )). To study B s — B s mixing and 
CP violation are very interesting topics and can be done only in hadronic environment. For 
experimental studies of B s —Bs mixing and certain CP violation in B s decays, the B s mesons 
being tagged at their production position are crucial. In hadronic environment to obtain 
copious B s mesons which are precisely tagged at their production positions, the way through 
B c decays is worth for considering seriously especially at LHC. In addition, potentially this 
way has a great advantage in rejecting backgrounds and deducing the systematic errors 
P). Therefore the feasibility for studying the B s mesons in this way should be investigated 
carefully, especially, by means of Monte Carlo simulating the events under specific detector 
environment so as to see the precise efficiency to have the tagged B S (B S ) mesons and whether 
the idea is really practicable or not. 

The paper is organized: following Introduction, in Section II we present the studies of 
the uncertainties in the estimates of the hadronic production for the dominant mechanism 
of gluon-gluon fusion and also a brief comparison between the mechanisms of gluon-gluon 
fusion and quark-antiquark annihilation. In Section III, we present the production with 
possible kinematic cuts, which roughly match the detector situation at TEVATRON and 
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FIG. 1: (a). The typical Feynman diagrams for the gluon-gluon fusion subprocess; (b). The typical 
Feynman diagrams for quark-antiquark annihilation subprocess. 

LHC. We also compute the production at TEVATRON with CM. energy changes from 
RUN-I to RUN-II. In the last section, we make some discussions and a short summary. 

II. THE UNCERTAINTIES IN ESTIMATES 

As stated in the introduction, there are two mechanisms for the lowest order of pQCD 
«f approach in the hadronic B c production, the gluon-gluon fusion and the quark-antiquark 
annihilation. Of them the gluon-gluon fusion mechanism is the dominant one, and there are 
36 Feynman diagrams for its subprocess. The typical Feynman diagrams for the subprocess 
are plotted in FIG. 1(a). In comparison, for the other mechanism, its subprocess is of quark- 
antiquark annihilation and there are 7 Feynman diagrams. The typical ones are plotted in 
FIG. 1(b). In fact, if we consider the production only through the dominant gluon-gluon 
fusion mechanism, then the contributions from the quark-antiquark annihilation mechanism 
will be considered as one kind of uncertainties. In this section, firstly we will focus on the 
gluon-gluon fusion mechanism and then we will make a brief comparison between these two 
mechanisms. 
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To be more useful experimentally, we study the uncertainties of B c hadronic production 
by means of the 'complete approach', i.e., at collision CM. energy \/S, 

da(S,p T ,---) = J " dx ± J dx 2 F l Hi Pi (x 1 , nl) ■ F 3 H2> p 2 (x 2 , fx 2 F ) 

ij 

■ da ij ^ BcbE (Pi,P 2 ,x 1 ,X2,iJ 2 F , s,pt,~-) > (!) 

where F % Hl P (xi, /j, f ) and F 3 H p 2 (x 2 , fJ-%) are the PDFs of incoming hadrons Hi (momentum 
Pi) and H 2 (momentum P 2 ) for parton i (with the momentum fraction xi) and parton j (with 
the momentum fraction x 2 ) respectively; \xf is the energy scale where the factorization and 
renormalization for the PDFs and the hard subprocess are made; dcfij-+B c be is the differential 
cross-section of the relevant hard subprocess, in which s = X\X 2 S is the c.m.s. energy of the 
subprocess and pt is the transverse momentum of B c . In the sub-sections below, we discuss 
the uncertainties which are from the non-perturbative factors and from those related to the 
hard subprocess. 



A. The Uncertainties Relevant to the Parameters of Potential Model and Masses 
of Quarks and the Meson 

For all the approaches and the color singlet mechanisms to the B c production, the 'decay 
constant' fs c relating the hadron matrix element, as one of important input-parameters, is 



needed as input, and so far it may be calculated either by potential model [1J| or by lattice 
QCD|l5|, while in the parer, we will make discussions related to it on potential model only. 
Within potential model, it directly relates to the wave function at the origin of the binding 
system in terms of the following formula, 



f , _ 12|¥ [1S] (0)| 2 



fl = ' 11J1V " , (2) 
m Bc 

where is the B c meson mass, \Pfis] is the wave function at the origin of the binding 
system (cb) at the [15] level. Since the spin splitting effects are ignored here, so there is no 
difference for the decay constant between the spin stats pSo] an d [ 3 Si]. There are several 
parameters in the potential model, which need to be fixed by fitting the experimental data 
of the heavy quarkonia ((cc) and (bb)), so when computing the spectrum and the wave 
functions of the (cb) svstem|l4j|. certain uncertainties from the fitting procedure will cause 
some uncertainties of /g c . In the leading order (LO) approximation for the production, f Bc 



TABLE I: The total cross section (in unit nb) for hadronic production of -B^l 1 ^] (B* [l 3 Si]) with 
various values of the c-quark mass m c and fixed 6-quark mass m& = 4.9 GeV. For definiteness and 
focussing the uncertainties from the quark masses alone, the rest uncertainty sources are fixed: the 
mass of the meson B c , m,B c = m c + m&; the gluon distribution function is taken from CTEQ5L; 
the factorization energy scale is chosen fj,p = Q 2 = s/4 and the running a s is of leading order. 





TEVATRON (VS = 2. TeV) 


LHC (v/5 = 14. TeV) 


m c (GeV) 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


a Bc ( nb ) 


4.84 


3.87 


3.12 


2.56 


2.12 


1.76 


75.6 


61.0 


49.8 


41.4 


34.7 


28.9 


a B * (nb) 


12.3 


9.53 


7.39 


5.92 


4.77 


3.87 


194. 


153. 


121. 


97.5 


80.0 


66.2 



appears in the amplitude as a linear factor exactly, so the production cross sections are 
proportional to it squared. Therefore, the uncertainties in the production from fs c can be 
figured out straightforwardly, so throughout the paper, we will fix the value /b c = 0.48 GeV 
for the ground states (the relevant uncertainties would be figured out if need). 

In addition to the decay constant fs c , the quark mass values m c and mj, also 'generate' 
uncertainties for the hadronic production. If all of the rest parameters in the potential might 
be fixed, then the values of quark masses would determine the decay constant /b c , the me- 
son mass mg c and etc. completely. However at present, the parameters in potential model, 
which also appear in the formula Eq.lfUl. cannot be completely fixed by fitting the available 
data of the heavy quarkonia, so the relations of the quark masses to the decay constant 
fs c , the meson mass and etc. cannot be well determined. Therefore, when studying the 
uncertainties for the hadronic production, we can consider all of the factors, explicitly ap- 
pearing in the potential model and in the formula for the subprocess Eq.([T|). in a 'factorized' 
way. Furthermore, since B c is the non-relativistic and weak-binding bound state, at LO the 
relative momentum between the constitute quarks can be ignored, e.g., we approximately 
have m,B c — rrib + m c . 

Throughout the paper, we study the uncertainties in 'a factorization way', i.e., all of the 
parameters vary independently in their reasonable regions. For instance, when focussing on 
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TABLE II: The total cross section (in unit nb) for hadronic production of -B c [l 1 5o] and £?*[l 3 Si] 
with various values of the 6-quark mass mj and fixed c-quark mass m c = 1.5 GeV. For definiteness 
and focussing the uncertainties from the quark masses alone, the gluon distribution function is 
taken from CTEQ5L; the factorization energy scale is chosen [ip = Q 2 = s/4 and the running a s 
is of leading order. 





TEVATRON (VS = 2. TeV) 


LHC (VS = 14. TeV) 


m b (GeV) 


4.5 


4.7 


4.9 


5.1 


5.3 


4.5 


4.7 


4.9 


5.1 


5.3 


<?B c ( nb ) 


4.10 


3.55 


3.12 


2.70 


2.38 


63.4 


56.2 


49.8 


44.1 


39.6 


a B * (nb) 


9.43 


8.35 


7.39 


6.59 


5.89 


148. 


133. 


121. 


110. 


100. 



the uncertainties from m c , we let it be a basic 'input' parameter varying in a possible range 

1.3GeV <m c < 1.8GeV, (3) 

with all the other factors, including the 6-quark mass, the decay constant fs c and etc. being 
fixed. 

The uncertainties from m c are indicated by the calculated total cross sections with the m c 
of Eq. Q , which are given in TABLE U Note that for the mass of _B C , the experimental result 
is itlb = 6.4±0.4 GeV [3], while the prediction by the potential model gives mg, = 6.1 ~ 6.3 
GeV Q and that by lattice QCD is about 6.4 GeV Q- Thus with m b = 4.9 GeV and 
rriB c — mf, + m c , the obtained m^ c is in the region of theoretical prediction as well as 
experimental measurement. The uncertainties from rrib can be analyzed in a similar way, 
where the 6-quark mass runs over the range: 

A.5GeV <m b < 5.3GeV. (4) 

The results for various b-quark mass are put in TABLE|TI] In Tables HI and HT1 the total cross- 
section for the hadronic production of ^[l 1 ^] and _B*[l 3 Si] at TEVATRON and LHC are 
computed, where the other factors are fixed precisely as: the PDFs are taking as CTEQ5L 
2(| ; the strong coupling a s is in leading order and the factorization energy scale is taken to 
be up = Q 2 = s/4, where s is the CM. energy squared of the subprocess. 

From the tables, one may observe how the values of m c and mb affect the cross section up 
to a sizable degree. Roughly speaking, both at TEVATRON and LHC, when nib increases by 



steps of 0.2 GeV, then the cross section decreases by 10% ~ 20%, while when m c increases 
by steps of 0.1 GeV, the cross section also decreases by 10% ~ 20%. 

In further investigations of the paper when examining the uncertainties from the rest 
factors, we shall take the values, m c — 1.5 GeV and = 4.9 GeV, for the quark masses. 



B. Uncertainties Relevant to the QCD Parameters and PDFs 

In this sub-section, let us focus on the gluon-gluon fusion mechanism only. 

As is shown in Eq.flTJ), PDFs F^^^xi, fi 2 F ) and F^ 2 , p Jx 2 , fJ,%) with if^proton, 
ff2=proton for LHC; iJi=proton, if2=anti-proton for TEVATRON generate uncertainties 
in the production. PDFs are of non-perturbative nature, and in Eq.Q they have been fac- 
torized out at the energy scale [ip with the help of pQCD factorization theorem. Usually 
the factorization is carried out at Q 2 , a characteristic energy scale for the hard subprocess, 
i.e. \xp = Q 2 , so we use Q 2 and /i F on the equal footing throughout the paper. 

As stated in the Introduction, all the estimates of the production are of leading (the 
lowest) order of pQCD and then how to choose the energy scale Q 2 is a very tricky problem. 
If Q 2 is chosen properly, the results may be quite accurate |31[. From experience, for a 
hard subprocess with a two-body final state, generally when taking Q 2 = |s ( is the 
CM. energy of the subprocess), we can achieve quite an accurate result for LO calculations. 
However, in the present case, the gluon-gluon fusion subprocess is of three-body final state 
and contain heavy quarks, so there are ambiguities in choosing Q 2 and various choices of Q 2 
would generate quite different results. Since such kind of ambiguity cannot be justified by 
the LO calculation itself, so we take it as the uncertainty of the LO calculation, although 
when the NLO calculation of the subprocess is available, the uncertainty will become under 
control a lot. While the NLO calculation is very complicated and it cannot be available 
in the foreseeable future, so here we take Q 2 as a possible characteristic momentum of the 
hard subprocess being squared. Having reviewed the choices of Q 2 in the literature, in the 
following we choose a few typical examples to study this kind of uncertainties. 

According to the factorization, the running of a s and PDFs should be of leading logarithm 
order (LLO), and the energy scale Q 2 appearing in the calculation should be taken as one 
of the possible characteristic ener gy s cales of the hard subprocess. In principle, as has been 
shown in references QflflHHHQ, there are several typical ways to choose the energy 
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scale Q 2 and to relate it to a characteristic energy scale (Q 2 3> A.q CD ) of the subprocess, 
such as Q 2 = 4m 2 , Q 2 = 4m^, Q 2 = Pr + rn 2 Bc (the 'transverse mass' squared of B c ), Q 2 = |s 
and etc.. When studying the uncertainties on the various choices of Q 2 , we precisely choose 
the following four types of Q 2 : 

Type A: Q 2 = s/4 , the CM. energy squared of the subprocess that is divided by 4; 

Type B: Q 2 = p\ + m 2 Bc , the transverse mass squared of the B c meson; 

Type C: Q 2 = s , the CM. energy squared of the subprocess; 

Type D, Q 2 = r pi rh + m\ , the transverse mass squared of the b quark. 

Since the PDFs can be obtained only through global fitting of the experimental data and 
evolute them to the requested characteristic scale Q 2 in standard way of pQCD, so there 
are several groups, CTEQ Q, Q, GRV H and MRS Q etc, who devote themselves 
to offer accurate PDFs to the world, and to keep PDFs updated as soon as possible when 
new relevant experimental data are available. Thus in literature, different versions of PDFs 
(including different issues by the same group) are used in the estimates of the hadronic 
production. To study the uncertainties due to various versions of PDFs, we consider the 
above mentioned three groups. 

The generator BCVEGPY is programmed based on the dominant subprocess to LO 
pQCD only, to be self-consistent, when applying BCVEGPY we shall adopt the LO PDF. 



Thus we take CTEQ6L 



28| and MRST2001L|29( as typical 



CTEQ5L [26], GRV98L 
examples for PDFs. The versions of the gluon distributions ended with 'L' are accurate 
up-to the leading logarithm order (LLO), i.e., their QCD evolution effects with a s running 
are included, so for the production to show the uncertainties correctly up to LO accuracy, it 
is necessary for the PDFs, the hard subprocess and the QCD 'coupling constant' a s 'run' to 
the energy scale Q 2 properly. The running a s in the PDFs and in the subprocess should be 
the same, and then when computing the production and taking the PDFs from one version 
of the three groups, the running a s should also be taken from the same group. 

For comparison between TEVATRON and LHC and to pinpoint the uncertainties from 
PDFs, a s running and the choices of the characteristic energy scale Q 2 , we calculate the 
production cross sections according to four types of PDFs, the strong coupling a s up to LLO 
fixed by the PDFs and the characteristic Q 2 chosen as Type A and Type B. The obtained 
results are shown in TABLE IIHI From the table, we can see that the gluon distribution 
functions, the running coupling constant and the corresponding energy scale, all affect the 
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TABLE III: Total cross-section for the hadronic production of -B^l 1 ^] and -B*[l 3 Si] at TEVA- 
TRON and at LHC with the leading order (LLO) running a s and the characteristic energy scale 
Q 2 = s/4 or Q 2 = pj, + m 2 Bc . The cross section is in unit of nb. 





CTEQ5L 


CTEQ6L 


GRV98L 


MRST2001L 


CTEQ5L 


CTEQ6L 


GRV98L 


MRST2001L 




Q 2 = s/4 


Q 2 =Pt + m 2 Bc 




TEVATRON 




3.12 


3.79 


3.27 


3.40 


4.39 


5.50 


4.54 


4.86 


a B*(l 3 S 1 ) 


7.39 


9.07 


7.88 


8.16 


10.7 


13.4 


11.1 


11.9 




LHC 


a B c (l^S ) 


49.8 


53.1 


53.9 


47.5 


65.3 


71.1 


70.0 


61.4 


°"fl«(l s Si) 


121. 


130. 


131. 


116. 


164. 


177. 


172. 


153. 



cross section to certain degrees. 

For all of the LO PDFs, the differences in the corresponding cross sections given by them 
are small (about 10% — 20%). But the choice of Q 2 as Type A and Type B causes the 
cross section to change by a factor of about | to |, which is a comparatively large effect. 
From TABLE llll| one may see that the cross section at TEVATRON and at LHC presents a 
somewhat different aspect of the cross section. Roughly speaking at TEVATRON CTEQ6L 
is the biggest, while at LHC, GRV98L is the biggest, under a certain characteristic energy 
scale such as Q 2 in type A. 

From TABLE IIII1 we may also see that the cross section of the B c meson production at 
LHC are at least one order larger in magnitude than that at TEVATRON. This is mainly 
due to the fact that LHC has much higher collide energy than TEVATRON, so the lowest 
boundary of the momentum fractions x« (i — 1,2) at LHC are much smaller than that at 
TEVATRON and then there are more interacting partons i.e. gluons, which have a CM. 
energy above the threshold for the subprocess, in the collision hadrons at LHC than at 
TEVATRON. 

When the produced B c mesons with a small pt are too close to the collision beam, 
they cannot be measured experimentally, so we calculate the cross-sections with transverse- 
momentum cuts pTcut- For comparison, we define the following ratio of the integrated cross 
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TABLE IV: Dependence of the ratio R = ( sxsxxbqil ) on p Tcut for B e [l 1 S Q ], where Q 2 = si A, the 
running a s is in LO and the gluon distribution is CTEQ5L. 





B c 


B* 


P T cut(GeV) 





5 


50 


100' 





5 


50 


100 


R(xlQ- 2 ) 


6.26 


5.24 


1.00 


0.31 


6.09 


5.19 


0.97 


0.24 




FIG. 2: B c differential distributions versus its transverse momentum px and rapidity y for different 
versions of gluon distributions in leading order. The characteristic energy scale is taken as Q 2 = s/4. 
Solid line: CTEQ5L; dotted line: GRV98L and dashed line: MRST2001L. The upper (lower) three 
lines corresponding to the distributions in LHC (TEVATRON). 

sections for TEVATRON and LHC, R = ( stmyatron ) anc i put the results for £ C [1%1 

V °LHC Jp Tcut * C L Uj 

in TABLE |IV| 

To show the uncertainties due to various factors more quantitatively, let us compute the 
differential cross sections of pt and y for the pseudo-scalar B c meson and draw the curves 
accordingly. 

From FIG. 121 we may see that the differential distributions for the three PDFs CTEQ5L, 
GRV98L and MRST2001L are very similar. The corresponding differences for TEVATRON 
and LHC are evident (the differences as shown in TABLE IlIII are less than 20%). In regions 
of comparatively small pt and y in TEVATRON the distributions of px and y obtained with 
MSRT2001L are the largest, and then those of GRV98L and CTEQ5L; while at LHC, those 
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FIG. 3: B c differential distributions versus its transverse momentum pr and rapidity y for four 
typical choices of the characteristic energy scale Q 2 . The gluon distribution is chosen as CTEQ5L 
and the running a s is in leading order. The choice of Q 2 is: solid line Type A; dotted line Type 
B; dashed line Type C and dash-dot line Type D. The upper (lower) four lines corresponding to 
the distributions at LHC (TEVATRON). 

of GRV98L are the largest, and then those of CTEQ5L and MSRT2001L. From the figure, 
we also see that the pt distributions in TEVATRON are steeper than those in LHC. 



as 



Since in 



10 



11 



iterature, various a s have been adopted in making estimates of the B c production 
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so to follow these authors, we try several a s to compute the production. 
Namely, in addition to the LO a s , the constant a s = 0.22 and the NLO a s are also taken into 
the computations and the curves of the production obtained by the various a s are drawn in 

fig. a 

We may see from FIG. El that the variations of the results are quite large for both TEVA- 
TRON and LHC; those with a s = 0.22 are the largest, those with NLO a s are the smallest, 
while those with LO a s are in the middle. 

In summary, of all the uncertainties which have been studied in this subsection, the 
uncertainty caused by the different choice of Q 2 is the largest, that can be as great as a 
factor around | (although in the literature results with various choices of a s running have 
given rise to even a greater change). 
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FIG. 4: B c differential distributions versus its transverse momentum px and rapidity y for different 
running a s . The gluon distribution of CTEQ5L and the Type A characteristic energy scale Q 2 = 
s/4, are used here. The solid line stands for the constant a s = 0.22, the dotted line for the leading 
order running a s and the dashed line for the next to leading order running a s . The upper (lower) 
four lines correspond to the distributions in LHC (TEVATRON). 

C. Comparison between the Mechanisms of Gluon-Gluon Fusion and Quark- 
Antiquark Annihilation 

In this subsection we make a brief comparison between the mechanisms of gluon-gluon 
fusion and quark-antiquark annihilation. Since, in a hadron such as P, P, ir, ■ ■ •, the PDFs 
of heavy quarks are very small throughout the active region \xi\ < 1, it is sufficient to take 
into account only the light quark pair (u-u, d-d and s-s) annihilations. 

The comparative results for the gluon-gluon fusion and the quark-antiquark annihilation 
mechanisms. In FIG we show 

From FIG we see that for production at a high energy, such as TEVATRON and 
LHC, the gluon-gluon fusion mechanism is dominant over the quark-antiquark annihilation 
mechanism, typically when the transverse momentum pt < 60 GeV and rapidity \y\ < 
5.0. Only in TEVATRON do the contributions from gluon-gluon fusion become smaller 
than those from quark-antiquark annihilation when pt > 70 GeV. Generally speaking, 
for the transverse momentum distributions, at low px the contributions from gluon-gluon 
fusion are greater than those from quark-antiquark annihilation by at least two orders of 
magnitude for both LHC and TEVATRON. For high transverse momentum (jp? > 70 GeV), 
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FIG. 5: B c distributions versus its transverse- momentum px and rapidity y of B c due to gluon- 
gluon fusion and quark-antiquark annihilation (the direction of positive y is defined as that of 
the momentum of the colliding proton). The PDFs are taken from CTEQ5L, a s running is set at 
leading order, and Q 2 = is. The solid lines correspond to gluon-gluon fusion and the dotted lines to 
quark-antiquark annihilation in TEVATRON. The dashed-dotted lines correspond to gluon-gluon 
fusion and the dashed lines to quark-antiquark annihilation in LHC. 

the difference in the two mechanisms changes sign in TEVATRON (i.e. contributions from 
quark-antiquark annihilation become greater than those from gluon-gluon fusion due to 
the fact that with such a high px, the valence quark PDFs play a significant role), but in 
LHC the contributions from gluon-gluon fusion are still greater by an order of magnitude 
than those from quark-antiquark annihilation. For rapidity y distributions, in the \y\ ~ 
region the contributions from gluon-gluon fusion are greater by two orders of magnitude in 
both colliders; but in TEVATRON in the higher region of 2.0 < \y\ < 5.0, contributions 
from the valance quarks' annihilations become more important, so that contributions from 
quark-antiquark annihilation increase remarkably (though still below those from gluon-gluon 
fusion). It is interesting to note that there is a slight asymmetry in the rapidity y of the 
meson B c for the quark-antiquark annihilation subprocess. In TEVATRON it is of PP 
collision, so the asymmetry in y is maintained somewhat (see FIG|3J) . We will describe the 
feature in TEVATRON elsewhere 0. 

The contributions from quark-antiquark annihilation mechanism to B c hadronic produc- 
tion, compared with those from gluon-gluon fusion mechanism, can be ignored safely in LHC 
in the allowed kinematic region. In TEVATRON, they can also be ignored over most of the 
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allowed kinematic region, i.e., only a little allowed kinematic region is an exception. In the 
paper, therefore, we consider the gluon-gluon fusion mechanism mainly. 

III. KINEMATIC CUTS FOR THE HADRONIC PRODUCTION OF B c MESON 

Experimentally, there is no detector which can cover all the kinematics of the events, so 
only some of B c produced events can be observed completely. For instance, in a high energy 
hadronic collider, the B c events with a small px and/or a large rapidity y (the produced 
B c mesons move very close to the beam direction) cannot be detected by the detectors at 
it directly, so this kind of events cannot be utilized for experimental studies in common 
cases. Therefore, only 'detectable' events should be taken into account in the estimates for a 
specific purpose, i.e., events with proper kinematic cuts on px and y must be put on precisely 
in the estimates. Considering detectors' abilities and to offer experimental references, we 
try various cuts accordingly in the estimate of the B c production. 

First, we study the distributions of px and \/I for the B c meson with various rapidity 
cuts y cut (here we mean that only the B c events with \y\ < y cut are taken into account), 
and as an example only the pseudo-scalar B c meson is taken into account. Note that the 
momenta and energies of the b, c and B c meson are measurable experimentally (a complete 
measurement of the events), hence the CM. energy \/§ of the subprocess is measurable, i.e., 
in fact the distribution of is measurable experimentally. Take into account the abilities in 
measuring rapidity of B c for the detectors CDF, DO and BTeV at TEVATRON, and ATLAS, 
CMS and LHC-B at LHC, various possible rapidity cuts, y cut ~ 1.5 or higher, are tried. To 
project out the cut effects, here we fix to use CTEQ5L for gluon distribution function, LO 
running a s and Type A energy scale (Q 2 = s 2 /4) to carry out the study. The results for 
the distributions of px and s with various y-cuts are shown a similar behavior, thus we 
do not plot the curve here, but note that the dependence of the differential distributions 
on y cut and pxcut for LHC is stronger than that for TEVATRON. This is because that at 
TEVATRON, all the distributions with sizable rapidity cover a smaller region in y. The 
correlations between px and y are interesting, so we plot the ^/-distributions with various 
Pr-cuts over a wide range pxcut '■ 5.0 ~ 100 GeV in FIG. El From the figure (FIG. IB}, the 
dependence of the differential distributions on rapidity y with different pxcut at LHC exhibits 
a broader profile than that at TEVATRON. 
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FIG. 6: B c differential distributions versus its y with various pTcut in TEVATRON (left diagram) 
and in LHC (right diagram). Solid line corresponds to the full production without PTcut] dashed 
line to PTcut = 5.0 GeV; dash-dot line to pxcut = 20.0 GeV; the dashed line to pxcut = 35.0 GeV; 
the big dotted line to PTcut = 50.0 GeV and the solid line with diamonds to PTcut = 100 GeV. 

The pT-distributions of the production vary with y cut mainly due to the fact that as 
p T increases, the dependence of the distribution on y becomes smaller as the value of y cut 
becomes less important. When p? increases to sufficient large value, the 'border' of the 
differential distributions (see FIG. \§§ becomes smaller than y cut and the pt differential 
distributions with and without y-cvX is coincide. 

To analyze the quantitative difference of the differential distributions with regard to prcut 
and y cut , we introduce a ratio for the integrated hadronic cross sections: 



where a Vcut and o"o are the hadronic cross section with and without y cut respectively. The 
ratio R PTcut varies with prcut and y cut , and its values are given in TABLE El 

From TABLE |V|we may see that, for a fixed y cut , the value of R PTcut becomes larger with 
increasing pxcut- It is understandable that the differential distributions versus the rapidity 
y decrease with the increment of p?, so the contributions to the hadronic cross section 
surviving after the cut, i.e. (\y\ < y cu t), increase with the increment of prcut- 

As for the B c meson production, only in high energy hadronic colliders can numerous 
B c mesons be produced. The cross section of B c meson production is quite sensitive to 
the CM. energy of the collider. Since RUN-II of TEVATRON is running now at a slightly 
higher energy than that 1.8 TeV (of RUN-I), we take TEVATRON as an example to study 




(5) 
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TABLE V: Values of the ratio R PTcut (see definition in text) for the hadronic production of pseudo- 
scalar B c meson in TEVATRON and LHC. 



PTcut 


0.0 GeV 


5 GeV 


20 GeV 


35 GeV 


50 GeV 


Dcut 


1.0 


1.5 


2.0 


1.0 


1.5 


2.0 


1.0 


1.5 


2.0 


1.0 


1.5 


2.0 


1.0 


1.5 


2.0 


Rp Tcut (TEVATRON) 


0.45 


0.64 


0.79 


0.46 


0.65 


0.80 


0.57 


0.77 


0.91 


0.65 


0.85 


0.95 


0.70 


0.90 


0.98 


R PTcut (LHC) 


0.31 


0.46 


0.59 


0.32 


0.47 


0.60 


0.38 


0.54 


0.69 


0.42 


0.60 


0.74 


0.45 


0.64 


0.79 



TABLE VI: The integrated hadronic cross section for TEVATRON at different CM. energies. The 
gluon distribution is chosen from CTEQ5L and the characteristic energy scale of the production 
is chosen as Type A, i.e. Q 2 = s/4. In addition, a cut for transverse momentum pt (pt < 5 GeV) 
and a cut for rapidity y (\y\ > 1.5) have been imposed. 



CM. energy 


1.8(TeV) 


1.9(TeV) 


1.96(TeV) 


2.0(TeV) 




0.40 


0.44 


0.46 


0.47 




1.00 


1.09 


1.14 


1.18 



the increase caused by slight changes in the CM. energy. The CM. energy for RUN-II is 
designed to be 2.0 TeV, so here we take energy values, 1.8, 1.9, 1.96 and 2.0 TeV to compute 
the cross sections. The obtained results for the values of the integrated cross-section with 
suitable cuts are shown in TABLE IV II From TABLE IVll we see that when the CM. energy 
of TEVATRON increases from 1.8 to 2.0 TeV, the cross-section increases by about 20%. 

IV. SUMMARY 

In this paper we have presented quantitative studies on the uncertainties in estimates of 
the B c meson hadronic production. The uncertainties have been examined in turn by 'fac- 
torizing' their origins. The computations are based on the generator BCVEGPY 23] mainly. 
The investigated quantitatively uncertainties involve those due to various versions of PDFs 
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given by various groups, the variations of the parameters relevant to the potential model, 
the strong coupling a s relevant to its running and the characteristic energy scale Q 2 of the 
process where the QCD factorization is carried out and etc.. We find that the characteristic 
energy scale Q 2 of the process is the source which causes the greatest uncertainty for the 
LO QCD estimates. We have also shown the differences between LHC and TEVATRON 
for various observables with reasonable kinematic cuts, such as the cuts on the B c meson 
transverse momentum pxcut and rapidity y cut . We also point out that at TEVATRON from 
RUN-I to RUN-II due to the increase of the collision CM. energy, the cross-section of the 
B c production increase by about 20%. In view of the fact that TEVATRON is running while 
LHC is still under constructing, based on our studies we can clearly see that experimental 
B c studies at the two colliders are stimulative and complimentary. As for the topic to study 
of B s — B s mixing and CP violations in B s meson decays at high energy hadronic collid- 
ers with tagged B s mesons produced through the B c decays, since LHC has much higher 
CM. energy, so the _B c -production cross-section is higher than one order of magnitude at 
TEVATRON, and LHC has much higher luminosity than at TEVATRON, it seems that the 
particularly interesting above topic may be more accessible and fruitful at LHC than that 
at TEVATRON RUN-II, even RUN-III. 
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